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ABSTRACT 

We investigate shattering and coagulation of dust grains in turbulent interstellar medium 
(ISM). The typical velocity of dust grain as a function of grain size has been calculated for var- 
ious ISM phases based on a theory of grain dynamics in compressible magnetohydrodynamic 
turbulence. In this paper, we develop a scheme of grain shattering and coagulation and apply 
it to turbulent ISM by using the grain velocities predicted by the above turbulence theory. 
Since large grains tend to acquire large velocity dispersions as shown by earlier studies, large 
grains tend to be shattered. Large shattering effects are indeed seen in warm ionized medium 
(WIM) within a few Myr for grains with radius a > 10~ 6 cm. We also show that shattering 
in warm neutral medium (WNM) can limit the largest grain size in ISM (a~2x 10~ 5 cm). 
On the other hand, coagulation tends to modify small grains since it only occurs when the 
grain velocity is small enough. Coagulation significantly modifies the grain size distribution 
in dense clouds (DC), where a large fraction of the grains with a < 10~ 6 cm coagulate in 10 
Myr. In fact, the correlation among the carbon bump strength, and the ultraviolet slope 
in the observed Milky Way extinction curves can be explained by the coagulation in DC. It is 
possible that the grain size distribution in the Milky Way is determined by a combination of all 
the above effects of shattering and coagulation. Considering that shattering and coagulation 
in turbulence are effective if dust-to-gas ratio is typically more than ~ 1/10 of the Galactic 
value, the regulation mechanism of grain size distribution should be different between metal- 
poor and metal-rich environments. 

Key words: dust, extinction — galaxies: ISM — ISM: evolution ISM: magnetic fields — 
methods: numerical — turbulence 



1 INTRODUCTION 

Dust grains absorb stellar ultraviolet (UV)-optical light and re- 
process it into far-infrared (FIR), the reby affecting the energet - 
ics of interstellar medium (ISM) (e.g. iHirashita & Ferrarall2002l) . 
For this process concerning the interaction between grains and ra- 
diation, the optical properties of grains are important. The grain 
optical properties are det ermined not only by dust species but 
also by the grain size (e.g. Dra ine & Ledl 19841) . In fact, the grain 
species and size distrib ution are derived f rom the observed inter- 
stellar extinction curve (Mathis e t alj 19771 . hereafter MRN). These 
grain properties also affect the FIR spectrum of dust emission (e.g. 
iTakeuchi et alj2005l) . 

The grain size distribution is known to be affected by 
various processes in the interstell ar space. Gr ains are sup- 
plied from stars at their death {G ehrz 198 ^) with a cer- 
tain grain size distribution toominik. Gail. & SedlmavJ 1 19891: 
Todini & Ferrara 2001; No zawa et alj|2003l ; iBianchi & Schneidea 
20071 : iNozawa et al.ll2007l) . These grains dispersed from stars are 



processed in the interstellar space. They grow in dense envi- 
ronments such as molecular clouds by grain-grain coagulation 
(Chokshi, Tielens, & Hollenbach 1993) and accretion of heavy el- 
ements IspitzeJ 1978b . They are also destroyed by supernova 
shocks by gas-grain sputtering and by grain-grain collision (shatter- 
ing) |Dwek & Scalo|l980l:lTielens et al.ll994l:lBorkowski & Dweld 
119951) . In particular. Uones. Tielens. & Hollenbach! d 19961 hereafter 
JTH96) show that the grain size distribution can be significantly 
modified by shattering. Such a change of grain size distribution 
would significantly modify the extinction curve and the infrared 
spectral energy distribution of dust emission. 

Potential importance of shattering and coagulation in ISM has 
often been pointed out, since relativ e velocity between grains is nat- 



urally expected if ISM is turbulent dKusaka et al.lll97C ;IV61k et al.1 
1 198a iDraindl 19851 : IOssenkopf|[~1993l : lLazarian & Yanll2002l). Be- 
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cause turbulence is ubiquitous in ISM (e.g. iMcKee & Ostrikeil 
l2007h , the relative grain motion induced by turbulence is of gen- 
eral importance in the grai n evolution in ISM . Moreover, the ISM 
is known to be magnetized jArons&MaxI 19751) . Thus, dust motion 
in magnetohydro dynamic (MHD) turbulence should be considered. 
lYan et al.l ( 120041 . hereafter YLD04) calculate the relative grain ve- 
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locity in compressible MHD turbulence, taking into account gas 
drag (hydr odrag) and gyro r esona nce. The basis of their theor y can 
be seen in|Lazarian & Yan (2002) and l Yan & Lazarianl d2003l) . Ac- 
cording to their results, grains can be accelerated to a velocity larger 
than a few km s _1 in diffuse medium by gyroresonance. At such a 
high velocity, grains can be shattered (JTH96). On the other hand, 
small grains can obtain velocities small enough for coagulation to 
occur, especially in dense medium (YLD04). 

The size distribution o f grains processed in vario us ISM 
phases was investigated by [GDonnell & Mafhjs| dl997l) . Their 
models incorporate coagulation by turbulent motion in clouds, ac- 
cretion of gas-phase metals onto grains, and shattering and sput- 
tering in interstellar shocks. They also take into account the phase 
exchange of the ISM. They show that both the observed extinc- 
tion curve and the observed depletion of refractory elements are 
reproduced by considering phase exchange among diffuse clouds, 
warm neutral intercloud gas, and molecular clouds. In their models, 
shatter ing mainly occurs in ISM shocks. However. I Yan & Lazarianl 
J2003h shows that grains can be accelerated to velocities large 
enough for shattering in MHD turbulence. Thus, it is worth fo- 
cusing on the effect of turbulence on the grain size distribution. 
In addition, the treatment of shattering can be revised by including 
the framework of JTH96 to take into account the velocity depen- 
dence of fragment production. It is also useful to compare the dif- 
ference between the size distribution modified by supernova shocks 
as treated by JTH96 and that by interstellar turbulence as examined 
in this paper. 

The aim of this paper is to examine quantitatively whether or 
not shattering and coagulation in turbulent ISM really modify the 
grain size di stribution. We focus on the effects of MHD turbulence 
as treated in lYan & Lazarianl fe003l) . The other types of dust pro- 
cessing such as shattering and sputtering in supernova shocks and 
dust condensation in stellar mass loss are not treated in this paper, 
in order to make our discussions concentrated and clear. For an ob- 
servation al comparison, we adopt the Milky Way extinction curve 
following lo'Donnell & Mathisl ll997l) . 

This paper is organized as follows. First, in Section[2] we de- 
scribe the model adopted to treat shattering and coagulation in tur- 
bulent ISM. Then, in Section [3] we overview the results. In Sec- 
tion [4] we discuss our results, focusing on the regulation mecha- 
nism of the grain size distribution in turbulent ISM. Section [5] is 
devoted to the summary. 



2 MODEL 

We consider the evolution of grain size distribution by shattering 
and coagulation induced by relative grain motions in turbulence. In 
our models, shattering and coagulation are treated simultaneously. 
The basic in gredients for shatteri ng and coagulation are taken from 
JTH96 and Cho kshi et al] (l993), respectively. We do not consider 
vaporization in grain collision, since this process does not change 
the grain size significantly at velocities (at most a few tens km s _1 ) 
achieved in turbulence. Our results are most sensitive to the col- 
lision rate between grains. Thus, first of all, the grain velocities 
adopted are discussed. Then, the frameworks of shattering and co- 
agulation are explained. 

2.1 Grain motion 

We assume spherical grains (Section [2.2| l. The velocity of a grain 
with radius a in the presence of interstellar MHD turbulence is 



taken from YLD04, who calculated the grain velocities achieved 
in various phases of ISM (CNM, WNM, WIM, MC, and DC, 
which stand for cold neutral medium, warm neutral medium, warm 
ionized medium, molecular cloud, dense cloud, respectively). The 
physical parameters which they adopted for each phase are listed in 
Table Q] For DC, they adopt two different cases for the ionization 
fraction (DC1 and DC2) to examine the effect of the uncertainty 
in the cosmic-ray ionization rate. We also adopt the same electrical 
charge of grains as calculated in YLD04, who considered photo- 
electric emission and collisions with ions and electrons based on 
the hydrogen number densities, the electron densities, and the UV 
radiation fields in Table[T] Since the values adopted for these quan- 
tities are those usually assumed for Galactic conditions, we expect 
that the calculations below are at least reasonable for the Galactic 
ISM. Grains are accelerated through turbulence hydrodrag and gy- 
roresonance. Below we briefly overview the models of turbulence 
and gyroresonance adopted by YLD04. 



2.1.1 Turbulence 

The velocity achieved by hydrodrag is determined by the largest 
scale on which grains are decoupled from the hydrodynamical mo- 
tion, since the turbulent velocity is larger on larger scales. Thus, 
larger grains, which tend to be coupled with larger-scale motions, 
are more accelerated. In general, neutral medium tends to acceler- 
ate grains less than ionized medium because of ion-neutral colli- 
sion damping of MHD turbulence. For example, the grain velocity 
in DC2 is smaller than DC 1 because of the difference in ionization 
degree. The same reason is applied to the larger grain velocities 
achieved in WIM than in WNM. 

YLD04 obtain the scaling relation of tur bulence velocity 
based on a MHD turbulence theory developed bv lCho & Lazarianl 
(2002), who decompose MHD fluctuations into Alfven, slow, and 
fast modes. Unlike hydrodynamic turbulence, Alfvenic turbulence 
is anisotropic, with eddie s elongated along the magnetic field 
( Gold reich & Sridhajl 1995h . This is because it is easier to mix the 
magnetic fields lines perpendicular to the direction of the mag- 
netic field rather than to bend them. The energies of eddies drops 
with the decrease of eddy size, and it becomes more difficult for 
smaller eddies to bend the magnetic field lines. Therefore, the ed- 
dies get more and more anisotropic as the sizes decreases. Eddies 
mix the magnetic field lines at the rate of k±v k , where k± is a 
wavenumber measured in the direction perpendicular to the local 
magnetic field and vt is the mixing velocity. The energy spec- 
trum for the perpendicular motions becomes Kolmogorov-like, i.e. 
Vk oc fc _1,/3 . On the other hand, the magnetic perturbations prop- 
agate along the magnetic field lines at the rate fcyVA, where fcii 
is the wavenumber parallel to the local magnetic field and Va is 
the Alfven velocity. The mixing motions couple to the wavelike 
motions parallel to magnetic field, giving a critical balance condi- 
tion, k±Vk ~ feii Va. Thus, we obtain fcii oc k 2 / 3 . The fast modes 
follow acous tic cascade, and show isotropic energy spectra with 
v k oc fc~ 1/4 dCho & Lazarianll2002l) . 

YLD04 assume that equal amounts of energy are transferred 
into fast and Alfven modes when driving is on large scales. 
The cascades proceed to small scales without much cross talk 
between those two kinds o f modes, according to the results in 
ICho & Lazarianl ( l2002l[2003l) . From the scaling relations of v k with 
k, we observe that the decoupling from fast modes usually brings 
larger velocity dispersions to grains than Alfven modes. 
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Table 1. The parameters of idealized ISM phases in YLD04. Among them, T is the gas temperature, rtjj is the 
number density of hydrogen atoms, n e is the number density of electrons, Guv is me UV intensity relative to the 
average local interstellar radiation field, L is the injection scale at which equipartition between magnetic and kinetic 
energies occurs, V is the effective injection velocity at the scale L (i.e. V = Va, where Va is the Alfven velocity), 
and k c is the wavenumber of the damping scale of turbulence. CNM: cold neutral medium; WNM: warm neutral 
medium; WIM: warm ionized medium; MC: molecular cloud; DC: dark cloud. 



ISM Phase 


CNM 


WNM 


WIM 


MC 


DC1 




DC2 


T (K) 


100 


6000 


8000 


25 




10 




njj (cm -3 ) 


30 


0.3 


0.1 


300 




10 4 




n c (cm -3 ) 


0.03 


0.03 


0.0991 


0.03 


0.01 




0.001 


Guv 


1 


1 


1 


0.1 


0.01 




0.001 


BQmG) 


6 


5.8 


3.35 


11 




80 




L(pc) 


0.64 


100 


100 


1 




1 




V = V A (kms" 1 ) 


2 


20 


20 


1.2 




1.5 




k c (cm -1 ) 


7 x 10~ 15 


4 X 10~ 17 




4.5 x 10" 14 


5.3 x 10~ 


15 5.3 


x 10~ 17 



In this paper, we do not consider imbalanced turbu- 
lence, which develops under unequal energy flux from oppo- 
site directions and ha s non- zero cross-helicity. Recent study by 
iBeresnvak & Lazarianl d2008l) shows that the stronger wave of the 
Alfven modes has smaller anisotropy, which indicates that the inter- 
action of the grains with the imbalanced Alfvenic turbulence could 
be more efficient. However, the results on the imbalanced turbu- 
lence is far from quantitative. And there is no conclusive theory yet 
for the imbalanced fast modes, which are more important for the ac- 
cel eration according to Yan & Lazarian ( 2003). The only study so 
far dSuzuki. Lazarian & Beresnvak||2007r) indicates that fast modes 
are not so different as in the balanced turbulence. In addition, im- 
balanced turbulence is applicable to places near an energy source, 
e.g. the vicinity of a star, which we do not consider in this paper. 

The conditions for ISM phases that YLD04 adopt imply that 
the turbulence is super-Alfvenic (SV > Va, where SV is the turbu- 
lence velocity). Indeed given the uncertainty of the strength of the 
magnetic field in the ISM, we do not know whether the conjecture 
is universal. However, our results for shattering and coagulation 
will not be concerned sensitively to the above debate on sub/super- 
Alfvenic turbulence. The reason is what follows. If the turbulence 
is sub-Alfvenic, the turbulence is weak. The weak turbulence has 
only limited inertial range. Moreover, it is the fast modes that dom- 
inate t he acceleration of dust as demonstrated by I Yan & LazariaJ 
d2003h . And there has been study showing that fast modes in weak 
turbulence are similar to fast modes in strong turbulence apart from 
the modes in the narrow cone around the k vector in Fourier space 
dChandrarJ[200l) . When the cascades proceed down to the scale 
where the critic al balance k±Vk ~ k\\VA is reached, turbulence 
becomes strong dLazarian & Vishniacll999i) . The coherence length 
of the strong turbulence LM A , is also the correlation length of the 
turbulence magnetic field, where L is the injection scale of the tur- 
bulence and Ma = SV/Va is the Alfvenic Mach number. It is un- 
likely that A/a is less than 0.1 in the Galactic environments. Given 
an injection scale of turbulence at 30 pc, the coherence length of 
the strong turbulence is then > 0.3 pc, which is still larger than the 
Larmor radius of most massive dust (a = 10 -4 cm). In fact, the 
observations ind icate that the correlation length of magnetic field 
is a few parsecs dSpangler & Mintonlll996T) . 

2.1.2 Gyroresonance 

Gyroresonance further accelerates charged grains. Grains obtain 
energy by resonant interactions with the waves if the resonance 



condition, u — knvfj, = nQ. (n = 0, ±1, ±2, ■ ■ ■), is satisfied, 
where ui is the wave frequency, ku is the parallel component of 
wavevector along the magnetic field, v is the particle velocity, fi 
is the cosine of the pitch angle relative to the magnetic field, and 
f2 = qB/(mc) is the Larmor frequency of the particle (q is the 
charge, B is the magnetic field strength, m is the grain mass, and 
c is the light speed). The above condition indicates that gyroreso- 
nance occurs when the Doppler- shifted frequency of the wave in 
the grain's guiding centre rest frame is a multiple of the gyrofre- 
quency and when the rotating direction of the electric wavevector 
is the same as the direction of the Larmor gyration of the grain. 
Then the steady state distribution function of grains is calculated 
by a Fokker-Planck equation treating the effects of gyroresonance 
acceleration and gaseous friction. 

Gyroresonance is efficient for large grains: The condition for 
gyroresonance is that the Larmor frequency Q. is smaller than the 
the cutoff wave frequency of the turbulence u c . The velocity of 
accelerated grain only weakly depends on the charge and mass as 
long as the aforementioned condition is satisfied (YLD04). On the 
other hand, the energy gain rate of the grains scales linearly with 
the intensity of MHD turbulence; thus the velocity is roughly pro- 
portional to the square root of the intensity of the MHD turbulence. 

2.1.3 Overall features of grain velocity 

The results in YLD04 indicate that gyroresonance accelerates sili- 
cate with a > 2 x 10~ 5 cm and graphite with a > 3 x 10 cm to 
velocities large enough (~ 20 km s _1 ) for shattering in WNM. 
Both silicate and graphite grains with a > several x 10~ 6 cm 
achieve velocities (1-2 km s _1 ) near to the shattering thresholds 
in CNM (Tabled. Although silicate and graphite with a > 10~ 6 
cm are accelerated to 20 km s -1 by gyroresonance in WIM, the 
acceleration by hydrodrag is larger because the dissipation of tur- 
bulence in WIM is less than that in WNM. 

On the other hand, low relative velocities of small grains al- 
low coagulation to occur. Moreover, with small velocities, a dense 
environment is necessary for a high enough collision rate. Thus, co- 
agulation is important in DC for grains with a < 10 -6 cm, which 
have velocities < 10 3 cm s~ 4 . Coagulation is also possible in MC 
with a smaller rate. 

Since YLD04 only calculated the grain velocity for a ^ 
10" 6 cm except for WIM, we extend the calculations down to 
a — 10 -7 cm. Below a ~ 10~ 6 cm, however, the coupling be- 
tween gas and grains occurs on a scale smaller than the dissipation 
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Table 2. Summary of grain properties. 



Species 


Pgr 


CO 


S «shat 


Pi 


Pv 


7 


E 


V 




(g cm~ 3 ) 


(kmr 1 ) 


(kms^ 1 ) 


(dyn cm -2 ) 


(dyn cm -2 ) 


(erg cm~ 2 ) 


(dyn cm -2 ) 




Silicate 


3.3 


5 


1.2 2.7 


3 X 10 11 


5.4 X 10 12 


25 


5.4 X 10 11 


0.17 


Graphite 


2.2 


1.8 


1.9 1.2 


4 X 10 10 


5.8 x 10 12 


12 


3.4 X 10 10 


0.5 



scale of turbulence. Thus, the velocities of grains typically smaller 
than 10~ 6 cm are determined by the thermal velocities. 



2.2 Time evolution of the grain size distribution 

We assume that grains are spherical with a constant material density 
Pgr. The mass m and the radius a of a grain are related by 



4tt 3 



(1) 



The number density of grains whose radii are between a and a + da 
is denoted as n(a) da, where the entire range of a is from a m i n to 
fflmax- The total grain mass is conserved in shattering and coag- 
ulation. To ensure the conservation of the total mass of grains, it 
is numerically convenient to consider the distribution function of 
grain mass instead of grain size. We denote the number density of 
grains whose masses are between m and m + dm as n(m) dm. 
The two distribution functions are related as n(a) da = n(m) dm. 

For numerical calculation, we consider TV discrete bins for the 
grain radius. The grain radius in the i-th (i — N) bin is 



between a' b \ and a\ a> , where a\°' — a^'-^S, a\°' = a m i n , and 

a^ = a max (i.e. log 5 specifies the width of a logarithmic bin: 
log 5 = (1/N) log(a max /a m i n )). We represent the grain radius 
and mass in the i-th bin with Oi = (afj 1 + a]° )/2 and m; = 
(47r/3)af Pg r . The boundary of the mass bin is defined as = 
(47r/3)[a' b ^] 3 p gr . Giving a m m, aimax, and TV, all bins can be set. 
A grain in the i-th bin is called "grain i". In this paper we take 
N — 32 after confirming that the results do not change if we take 
a larger TV. For the size range, we assume a m in = 0.001 pm and 
cimax = 0.25 pm to reproduce the Milky Way extinction curve 
(Sectionl2~4ll. 

The mass density of grains contained in the i-th bin, pi, is 
defined as 



(b) 



,00 



,0>) 



Pi = mih(mi)(mf^ - m^) . 
Then, the time evolution of pi is expressed as 



dpi 
dt 



dpi 
dt 



+ 



dp, 



dt 



(2) 



(3) 



where the first and the second terms in the right-hand side are the 
contributions from shattering and coagulation, respectively. These 
two terms are estimated in Sections [2.2.1l and l2.2.2l 

We consider silicate and graphite as grain species. In order 
to avoid complexity in compound species, we only treat collisions 
between the same species. Although this underestimates the grain 
collision rate by a factor of ~ 2, our simple assumption here is 
enough to understand the effects of interstellar turbulence on the 
grain size distribution for the first time. The adopted parameters for 
each grain species are summari zed in Tableland are taken from 
JTH96 and Cho kshi et al] |l993). We use the same notation (n(a)) 
for both silicate and graphite size distributions. 



2.2.1 Shattering 

The time evolution of pi by shattering can be written as 

JV WW 



dt 



liPi /J CtkiPk + /J / ] 



a kj pkPjm sh!it (i) , 



j = l k = l 



ifv ki > Ujhat, 



Ctki 



OhiVhi 



otherwise, 



(4) 



(5) 



where m*jj t (i) is the total mass of the shattered fragments of a 
grain ft that enter the i-th bin in the collision between grains k and 
j, Oki and Vki are, respectively, the grain-grain collisional cross 
section and the relative collision speed between grains k and i, and 
fshat i s tne velocity threshold for shattering to occur. For the cross 
section, we apply ati = ir(ak + at) 2 . 

The grain velocities given by YLD04 are typical velocity dis- 
persions. The relative velocity vm is treated with a similar manner 
to Appendix A of JTH96. Each time step is divided into 4 small 
steps, and we apply Vik = Vi + v^, \vi ~ vu\,Vi, and Vk in each 
small step, where Vi and Vk are the velocities of grains i and k, 
respectively (see Section [27Tb . Note that the mass distribution of 
the shattered fragment m^ at (i) depends on Vkj. The method for 
kj 



calculating m sl ^ at (i) is described in Section l23 



2.2.2 Coagulation 

The time evolution of pi by coagulation can be written in a similar 
form to equation (01 : 



dp; 



dt 



MPi ^2 a ^iPk + ^ ^ Ctfcj/5fcPj"l?oag(2) • 



j = l k = l 



a l: , 



OkiVki 

miirik 



if V ki < «coag, 

otherwise. 



(6) 



(7) 



Here, m^L 



, k 'J 



= mi if m£l\ ^ mfc + rrij < m' b '; otherwise 
(i) = 0. The coagulation is assumed to occur only if the rel- 
ative velocity is less than the coagulation thre shold velocity «^ ag . 
The coagulation threshold velocity is given by dChokshi et al J 19931 : 
iDominik & Tielensj|l99l YLD04) 



= 2.14F stic 



(ah 



1/2 



v 5/6 



E^R^p 



5/6 1/2 
gr 



(8) 



where a factor i^tick = 10 is introduced following YLD04 (based 
on the experimental work bv lBlurrj l2000). 7 is the surface energy 
per unit area, Rki = akai/(ak + Oj) is the reduced radius of the 
grains, E is related to Poisson's ratios (i>h an d fi) and Young's 
modulus (E k and Ei) by 1/E = (1 - VkfjEk + (1 - Vif/Ei. 
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The v alues of 7 and E are taken from Table 3 of IChokshi etal] 
Jl993h (the data for quartz and graphite are used for silicate and 
graphite, respectively) as summarized in Tabled The treatment of 
Vki in coagulation is the same as that in shattering (Section|ZZT). 

Although the above form of coagulation threshold velocity is 
derived based on both physical and experimental basis, there could 
be significant uncertainties. A change of the coagulation threshold 
affects the largest size of grains subject to coagulation, since larger 
grains have larger velocities. For example, as shown later, the coag- 
ulation condition is satisfied for grains with a < a few x 10 -6 cm. 
However, we have confirmed that even if the coagulation thresh- 
old is altered by a factor of ten from the above values, the size of 
grains subject to coagulation changes only by a factor of 2-3. This 
is because of a steep dependence of the grain velocity on the grain 






1 grain j (projectile) 


grain k (target) 


ify/t/>u stal 


1 if m k > Mil 


ifm k <M/2 | 








grain/ 



grain k 

mfrag = A4j m ftag = m > 



IN*? 



grain k 

'"fag = "« Wfag = mj 



2.3 Production of shattered fragments 

Here we determine the mass distribution of the shattered fragments, 
m shat(*)> m equation {4j». The shattering rate is determined mostly 
by the collision frequency between grains. Thus, the overall results 
in this paper is less sensitive to the detailed model of shattered frag- 
ments than to the grain velocities. Indeed, JTH96 show that the 
size distribution of shattered fragments does not have a large influ- 
ence on the overall grain size distribution. We have also confirmed 
that the specific parameters adopted in this subsection do not affect 
the results significantly. We can have an idea about the uncertain- 
ties caused by the material parameters by comparing the results for 
silicate and graphite. Since the results in these two materials are 
broadly similarly the assumption on the material parameters only 
has a secondary importance in our models as long as there are not 
very eccentric materials involved. 

We illustrate our treatment of shattering in Fig.Q] We consider 
a collision between grains k and j (here we assume k j), and call 
the grains labeled as k and j target and projectile, respectively. The 
necessary material quantities are summarized in Table[2] The mass 
shocked to the critical press ure for cratering in the target, M, is 
given by dTielens et all 1994 JTH96f| 



M 



1 + 211 



1 



(1 + CT l/9 



M r 



8/9 



(9) 



where 1Z = 1 in the collision between the same species (we adopt 
1Z = 1 in this paper), Mr = vtj /en (co is the sound speed of the 
grain material), <ti and an are constants typically of order unity 
(equation II lb . and Mi is the Mach number corresponding to the 
critical pressure Pi : 



Mi 



20i 



l + (l-r-4s</>i) 1 /2 



(10) 



where <j)i = Pi/(p gr c ), and s is a dimensionless material con- 
stant that determines the relation between the shock velocity and 
the velocity of the shocked matter. Using the following expression 
for a as 



a(M) 



0.30(s + M~ L - 0.11)° 



(11) 



s + M- 1 -! 

we obtain ax = a^Mi) and an = a(M r /(l + TV)). We assume 



Figure 1. Schematic figure of our treatment of shattering. We consider a 
collision of two grains in the fe-th size bin and the j'-th size bin (called 
grain k and grain j, respectively) with a relative velocity of ujy, We call 
the larger and smaller grains "target" and "projectile", respectively. Here 
we assume k j; that is, the target is grain k and the projectile is grain j. 
If the shocked mass M in the target is larger than the mass of the target 
(M > rrifc/2), we assume that the entire target fragments into small pieces 
( m frag = f^k)' If M < mfe/2, apart of the target mass, M e j is shattered 
and ejected (mf rag = Af c j). The entire projectile is assumed to fragment 
( m frag = m j )- The size distribution of the fragments is assumed to follow 
equation H2\ . 



that if more than the half of the target is shocked (i.e. M > rrik/2) 
the entire target is shattered (i.e. mf rag = rriu in equation [13] 
where mf rag is the total mass of the fragments). Otherwise, only 
a fraction of the target mass (M e j) is ejected from the target (i.e. 
m fra g = Af e j). M e j is assumed to be 0.40M, i.e. 40% of the 
shocked mass is finally ejected. This fraction is derived for z = 3.4, 
where the radial velocity of the cratering flow in the shattered ma- 
terial is approximated to be cx R~ z (R is the distance from the 
cratering centre; ITH96). Finally, the entire projectile is assumed 
to fragment into small pieces (i.e. mf rag = rrij for all projectiles). 

The frag ments are assumed to follow the size distribution 
faellveJl97d ; JTH96) 



"fr; 



, a 



-3.3 



da , 



trag(a) da — Cfra 
where the normalization constant Cf rl 



is determined by 



m fri 



?™ifrag(a) da . 



(12) 



(13) 



Here, af m i n and af max , respectively, specify the upper and lower 
bounds of the fragment radius, which are determined in Sections 
12.3.11 and 12.3.21 for the projectile and the target, respectively. If 
afmin estimated is less than a m in, we take afmin = a m i n . If af max is 
also smaller than a m i n , all the fragments are put in the bin with the 
smallest size (i.e. i = 0). Finally, the mass of shattered fragments 
in the i-th bin is determined in terms of nf rag as 



m s h a t(«) = / mn frag (a)da, 

„< b > 



(14) 



1 If the results are very different between silicate and graphite, the differ- 
ence can be attributed to the difference in grain velocity. 

2 In ITH96, ay is denoted as an. 



where the integration is performed in the range corresponding to the 
i-th bin. We put this mass in equation l[4j (note that the superscript 
kj is omitted here). 
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In the following, we summarize how to determine af m i n and 



2.3.1 Projectile rrij 

The entire projectile is assumed to fragment into small pieces, i.e. 
wifrag = rrij. The maximum grain size, af max , is determined by 
(JTH96) 



= 0.22a, 



"cat 



where u ca t is the critical spalling collision velocity given by 

9/16 



«cat = Co 



1/2 1/16/, . t-,s K , 



(15) 



(16) 



(1 + 2TZ)mj 
The minimum grain size, af m i n , is determined by 

afmin = 0.03af max . (17) 

2.3.2 Target m k 

If M/rrik > 0.5, we assume that the entire target fragments into 
small pieces; i.e. mf rag = m^. The maximum and minimum frag- 
ment sizes are determined by equations l |15l > and (TTJl, respectively, 
but j and k are exchanged. 

If M/m k < 0.5, we assume that the mass M ei (= 0.40M) 
fragments into smaller grains (i.e. raf rag = M e j), and a grain with 
mass of mfe — M e j is left, which is put in the corresponding bin. Ac- 
cording to equations (10) and (15) in JTH96, the largest fragment 
size and the total ejected volume (A/ c j/p gr ) can be related by 



Mei/p g 



16 z\z-1 



' 3 

" ^fmax ■ 



(18) 



We determine a mlax according to this equation with z — 3.4. The 
following estimate is adopted for af m i n (JTH96): 

, Pl xl.47 



m ■ 

where P v is the critical pressure for vaporization (TableO. 



2.4 Initial Grain Size Distribution 

It is not an easy task to select a good initial grain size distribution, 
since the grain production in stellar mass loss is not fully under- 
stood yet. Thus, we concentrate on how the standard grain size dis- 
tribution is modified by shattering and coagulation in various ISM 
phases. As the standard grain size distribution, we adopt 



n(a) = Ca (a min < a < a max ) . 



(20) 



Li & Draine 



where C is the normalizing constant. We select K = 3.5 as derived 
by MRN to explain the observed Milky Way extinction curve. For 
the size range, we assume a m i n = 0. 001 /tm and a max = 0.25 /J,m 
for both graphite a nd silicate (MRN: IlT~& Draine|[200lh . although 
^IfeOOlh adopt more elaborate functional form (see also 
in. & H endi -vl 19941) . In fact, as shown later, the predicted 



extinction curve is broadly consistent with the observed extinction 
curve (Section r3.2t . Thus, the above simple assumption for the size 
distribution is enough for our purpose. 

The normalization factor C is determined according to the 
mass density of the grains in the ISM: 



4?r 3 r -K A 
—a p gr Ca da , 



(21) 



where tih is the hydrogen number density given for each ISM 
phase in YLD04 (see also Table QJ, mn is the hydrogen atom 
mass, and 1Z is the dust-to-hydrogen mass ratio (i.e. dust abun- 
dance relative to hydrogen) in the ISM. We adopt 1Z = 4.0 x 
10~ 3 and 3.4 x 10~ 3 for silicate and graphite, respectively 



( Takag iTVansevicius. & Arimotoll2003h . As shown in Section [3~2l 
the size distribution assumed here reproduces the observed Milky 
Way extinction curve. 



2.5 Timescales 



We calculate the change of grain size distribution in various ISM 
phases on typical timescales. The typical timescale of the phase 
change among WIM, WNM, and CNM is ~ a few x 10 7 -10 8 
vrjlkeuchill 19881 ; lO'Donnell & Mathis|[l997l ; Iffirashita & Kamaval 
l200lE 7 For WIM, there is another relevant timescale, that is, recom- 
bination timescale. With hydrogen number density ~ 0.1 cm~ 3 
and tempe rature ~ 10 4 K, the recombination timescale is roughly 
~ 10 6 yr (Spitzer] [l978l) . As shown later, the grains are shattered 
too much in WIM for a time t > 10 7 yr, so a timescale of the 
order of Myr is more appropriate for WIM (Section[3j. For denser 
medium, a short lifetime may be reasonable, and indee d the lifetime 
of molecular clouds is estimated to be ~ 10 7 yr ([Blitz & Sh vjl98d; 
IPalla & Stahler 2002; Kaw amura et alj2007h or shorter ( Elmegreen 
2000; Hartmann 2003). Thus, we examine t < 10 7 yr for MC and 
DC. T hese timescales are also consistent with lO'Donnell &^ Mathis 
Jl997l) . 



2.6 Extinction curves 



Following lO'Donnell & Mafhisl Jl997h . we use extinction curves 
to test our results. We calculate extinction curves by using the 
optical constants of astronomical silicate and graphite taken from 
iDraine & Led d 19841) . Then cross sectio ns for absorption and scat- 
tering are calculated with Mie theory (Bohren & Huf fmanl [19831) 
and weighted for the grain size distribution. Finally the extinction 
curves of silicate and graphite are summed up. The extinction is 
normalized to the number density of hydrogen atoms. 

For comparison, the observational data of the standard in - 
terstellar extinction of th e Milk y Way is taken from IPeil i 19921) . 
iBohlin. Savage, & Drake] ( 1 19781) show that the mean Milky Way 
Nh/E(B — V) (Nh is the column density of hydrogen atoms 
and E(B - V) is the excess of B - V colour) is 5.8 x 10 21 
atoms cm~ 2 mag -1 . Then by using Ab = (1 + Rv)E(B — V) 
{A\ is the extinction in units of magnitude at wavele ngth A an d 
Rv = A V /E{B ~ V)), and adopting R v = 3.08 Jpellll992l) . 
we obtain Nn/A B = 1.422 x 1 21 a toms cm" 2 mag -1 for the 
mean Milky Way extinction. |PeH dl992h lists f (A) = A\/A B for 
relevant wavelengths, and the equation A\/Nh — £(\)Ab/Nh 
can be used to obtain A\/Nh for the mean Milky Way extinction 
curve. The extinction curves are often normalized to the value at V 
band, but we do not adopt this normalization, because the V band 
extinctions themselves in our models are significantly affected by a 
slight change of the size distribution at 10~ 6 cm < a < 10 -5 cm 
(Section[3]l. Since our models are based on a simple analytical treat- 
ment of interstellar turbulence with a single density, it is not reason- 
able to adopt a normalization parameter which is not robust to the 
change of details in the models. 

It is also known that there is a variation in the 
Milky Way extinction cu r ves a long various lines of sight. 
ICardelli, Clayton. & Mathisl d 19891) argue that the variation of the 
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extinction curves can be parametrized by Ry- More recently, 
iFitzpatrick & Massal fe007l) show that the variance of the extinc- 
tion curves normalized to Ay is roughly 20% at 1/ A = 8 fxm~ 
and roughly 10% at the 2175A bump. Although we do not adopt 
the normalization at V band in the extinction curve, these variances 
provide us with a rough idea as to how much variation of the ex- 
tinction curve is permitted in the Galactic environment. 



3 RESULTS 

3.1 Grain size distribution 

In Figs. [2] and [3] we show the size distributions of silicate and 
graphite, respectively, for various ISM phases. The size distribution 
is expressed by multiplying a 4 to sho w the "mass distribution" i n 
each logarithmic bin of the grain size JO'Donnell & Mathisll 19971) : 
i.e. a 3 comes from the grain mass and another factor a originates 
from da/d In a = a. The largest change is seen in WIM, for which 
we present the grain size distributions in the shortest timescales 
(1 Myr and 5 Myr). If the grains are processed for a longer time 
in WIM, the extinction curves become too modified to be consis- 
tent with the observed Milky Way extinction curve (Section [3.21 ). 
In WIM, grains with a > a few x 10~ 6 cm are efficiently acceler- 
ated by up to a velocity larger than the shattering threshold veloc- 
ity. If the grain velocity is the same, shattering efficiently destroys 
small grains because of their large surface-to-volume ratios. Thus, 
the largest shattering efficiency is realized for the smallest grains 
which obtain a velocity above the shattering threshold. For this rea- 
son, grains with a ~ 10 -6 cm are the most efficiently shattered in 
WIM. This is different from shattering in supernova shocks, where 
such small grains are not efficiently shattered (JTH96) since small 
grains tend to be decelerated quickly by the gas drag. 

In WNM, because of the ion-neutral collision, the fast modes 
are damped on a larger scale than in WIM. Thus, gyroresonance 
is not efficient for small grains, and only large grains with a > 
0.2 /im for silicate and with a > 0.3 /im for graphite can be accel- 
erated to a velocity large enough for shattering. (Note that graphite 
grains are not shattered since we only consider a ^ 0.25 fim.) 
Those threshold radii for gyroresonance are quite robust because 
they only weakly depend on the charge, the magnetic field strength, 
and the grain density as mentioned in Section l2.1.2l It is interesting 
that those grain sizes satisfying the shattering condition in WNM 
are nearly the upper grain size in the Milky Way (MRN). Thus, the 
upper limit of the grain size is possibly determined by shattering in 
WNM. This issue is further investigated in Section POl 

Shattering takes place also in CNM for graphite because 
graphite has lower shattering threshold velocity than silicate. How- 
ever, the result is sensitive to slight changes in the shattering thresh- 
old. Moreover, the grain velocities acquired by the gyroresonance 
have uncertainties coming from the magnetic field strength and the 
grain charge, although the uncertainties are generally small (Sec- 
tion [27TJ. Thus, the shattering in CNM is not conclusive. We ob- 
serve slight coagulation of silicate with a < 10 -6 cm in CNM. 

In MC and DC, coagulation takes place. In particular, an ap- 
preciable amount of small grains coagulate in DC because of high 
density. Since the grain velocities are lower than the coagulation 
threshold for a < a few x 10 -6 cm, the grains accumulates around 
a ~ a few x 10~ 6 cm in DC. Coagulation occurs up to larger grain 
radii in DC2 than in DC1 because the velocity is lower in DC2 than 
in DC1 because of ion-neutral damping of turbulence (Section ETTb . 



3.2 Extinction Curves 

For observational comparisons, we calculate the extinction curves 
with the method described in Section |2~6l In WNM and MC, the 
grain size distributions are modified too slightly to change the ex- 
tinction curves significantly. The interesting cases are WIM, CNM, 
and DC, for which we show the results in Fig. [4] First of all, the ini- 
tial MRN distribution reproduces the Milky Way extinction curve 
including the UV slope and the 2175 A bump. The only devia tion 
is seen at 1/A ~ 6 /im~ . The same deviation is also seen in 
(1992). Since the aim of this paper is not precise fitting of the ex- 
tinction curve, we do not fine-tune the grain size distribution. Ex- 
amples of detailed fit ting of the extinction curve can be seen in 
iKim et all Jl994h and lWeingartner & Drainel feOOll) . Below we de- 
scribe some features in the extinction curves calculated for WIM, 
CNM, and DC. 

3.2.1 WIM 

The extinction curves of the grains processed in WIM are shown 
in Fig. [4^. At t = 5 Myr, the 2175 A bump is too high and the 
UV slope is too steep to be consistent with the observed Milky 
Way extinction curve. Thus, we can conclude that the grains are 
continuously processed in WIM for no longer than 5 Myr. It is in- 
teresting to point out that this timescale is roughly comparable to 
the recombination timescale of g as with density ~ 0.1 cm -3 and 
temperature ~ 10 4 K (~ 10 B vr; ISpitzedri978n and to the typical 
lifetime of massive stars (source of ionizing photons). 

3.2.2 CNM 

In Fig. |4p, we show the extinction curves in CNM. Because large 
graphite grains are shattered, the 2175 A bump becomes high and 
the UV slope becomes steep. The extinction curves after shatter- 
ing in CNM do not deviate largely from the observed Milky Way 
extinction curve within 10 Myr. If grains suffer a longer time of 
shattering, the extinction curve shows too high a 2175 A bump and 
too steep a UV slope to be consistent with the observed extinc- 
tion curve. However, as noted in Section [37X1 the arguments here 
are sensitive to the prediction of grain velocities and the assumed 
value of shattering threshold velocities. 

3.2.3 DC 

In Figs. [4J; and d, we present the extinction curves in DC1 and 
DC2, respectively. Because the grain size is biased toward large 
sizes after coagulation, the 2175 A bump is lower and the UV slope 
is less steep than the original curve predicted from the MRN size 
distribution. 

The extinction curves in DC2 are more consistent with the 
observed Milky Way extinction curves than those in DC1 in the 
following two points. First, the wavelength at the peak of 2175 
A bump changes less in DC2 than in DC1. The observed central 
wavelengths of 2175 A bump in various lines of sig ht in the Milky 
Way are insensitive to the variation of bump strength iCard elli et al.l 
Il989l) . The different behaviours of the 2175 A bump between DC1 
and DC2 come from the "smoothness" of the size distribution 
around a ~ a few x 10~ 6 cm: In DC1 the graphite size distri- 
bution show a very steep depletion of grains at a < 3 x 10~° cm, 
while in DC2 the depletion of such small grains is not so drastic as 
inDCl. 

Second, the behaviours of Rv in terms of the 2175 A bump 
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Figure 2. The grain size distributions for silicate in (a) WIM, (b) WNM, (c) CNM, (d) MC, (e) DC1, and (f) DC2. The initial MRN distribution is shown by 
the dotted line. The solid (dashed) line presents the distribution at t = 1 Myr (t = 5 Myr) for Panel (a), t = 10 Myr and 50 Myr for Panels (b) and (c), and at 
t = 5 Myr (t = 10 Myr) for Panels (d), (e), and (f). In this paper, the grain size distributions are presented by multiplying a 4 to show the mass distribution in 
each logarithmic bin of the grain radius. The grain size distribution is normalized to the hydrogen number density ripj. 
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Figure 3. Same as Fig.|2]but for graphite. 
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Figure 4. The extinction curves of grains processed in various ISM phases. In each panel, the dotted line represents the initial MRN distribution, (a) t = 1 Myr 
and 5 Myr in WIM; (b) t = 10 Myr and 50 Myr in CNM; (c) t = 5 Myr and 10 Myr in D C1; an d t = 5 Myr and 10 Myr in DC2 for the solid and dashed 
lines, respectively. The squares show the data of the observed Milky extinction curve bv lPeil d 19921) . In this paper, the extinction is normalized to the hydrogen 
column density. 




Figure 5. The grain size distributions of (a) silicate and (b) graphite after t = 3 Myr (dashed line) and t = 5 Myr (dot-dashed line) in WIM from the initial 
distribution (solid line), for which the size distribution processed for 10 Myr in DC2 is adopted. The dotted line in each panel shows the MRN distribution. 



Dust shattering and coagulation in interstellar turbulence 



and the UV slope are more consistent with the observed Milky Way 
extinction curves in DC2 than in DC1. The observed Milky Way 
extinction curves show that a lar ge Ry is related to a weak 2175 A 
bump and a shallow UV slope JCardelli et al.ll 19891) . Starting from 
3.6 at t = 0, Rv changes to 3.2 (t = 5 Myr) and 3.1 (t = 10 Myr) 
in DC1, while it changes to 4.2 and 4.8 in DC2. Thus, DC1 has a 
trend opposite to the observed one, while DC2 reproduces a right 
trend between Ry, the 2175 A bump, and the UV slope. Only the 
grains with a < 4 x 1CP 6 cm, whose velocity is below the coagula- 
tion threshold velocity, can coagulate in DC1. Since the grain size 
change in this small size range does not affect the extinction in long 
wavelengths such as B and V bands0 coagulation of larger grains 
is necessary to change Rv - This is why Ry changes only a little in 
DC1. In DC2, coagulation to larger grain sizes indeed occurs and 
Ry increases as coagulation proceeds. Considering that there are 
uncertainties in the threshold velocity for coagulation and in the 
grain velocities, the success in reproducing qualitatively the trend 
among Ry, the 2175 A bump, and the UV slope in DC2 supports 
the view that coagulation induced by turbulent motions in dense 
environments really occurs in the ISM. 



4 DISCUSSION 

The important features found in the previous section can be sum- 
marized as follows. 

(i) The largest effect of shattering is seen in WIM, where 
grains with a > 10 -6 cm are efficiently shattered. 

(ii) The largest effect of coagulation is observed in DC around 
a < a few x 10 -6 cm. 

(iii) Grains with a > a few x 1CP J cm can be shattered in 
WNM and graphite grains with a > 10 cm may be quite effi- 
ciently destroyed in CNM. These destructions could affect the up- 
per limit of grain size in ISM. 

(iv) On the other hand, the lower limit of grain size may be 
determined by coagulation in DC and MC. 

The features (i) and (ii) indicate that once grains are included 
in WIM or DC, the grain size distribution is significantly modified. 
It is interesting to note that the shattered grains in (i) and the coagu- 
lated grains in (ii) have a similar size. Thus, it is worth investigating 
if the MRN size distribution can be realized as a balance between 
(i) and (ii). This point is investigated in Section |4~T| 

Regarding the feature (iii), as mentioned in Section 13.2.21 
the results in CNM are sensitive to the grain velocities and the 
shattering thresholds. We leave more careful treatment of shatter- 
ing in CNM for future work. Shattering in WNM is interesting 
to investigate, since turbulence in WNM accelerates grains with 
a > a few x 10~ 5 cm much above the threshold velocity for shat- 
tering. This size really matches the upper limit of the grain size 
distribution (MRN). This point is investigated in Section l4~2l 

The issue (iv) has already been investigated and discussed in 
Section EH^I 

4.1 Grain size distributions in diffuse-dense phase exchange 

In ISM, mass is ex changed betw een various phases 

dMcKee & Ostrikej 1 19771 ; llkeuchil 1 19881) . Thus, it is impor- 

3 From the knowledge of Mie theory, if the grain size is much smaller 
than X/2n, the extinction becomes inefficient, i.e. Q < 1, where Q is 
the extinction cross sectio n normalized to the geometrical cross section 
iBohren & Huffmanll 19831) . 
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Figure 6. The extinction curves calculated for the size distributions in Fig. 
\5\ The solid line represents the initial size distribution (10 Myr in DC2), and 
the dashed and dot-dashed lines show the extinction curves at t = 3 Myr 
and 5 Myr in WIM, respectively. 



tant to investigate the effects of multi-phase ISM on the evolution 
of grain size distribution, although the main aim of this paper is 
to examine the dust processing in individual phases. The largest 
shattering and coagulation effects are seen in WIM and DC, 
respectively, and we here examine the dust processing in both 
WIM and DC to address a possible importance of multi-phase ISM 
in determining the grain size distribution. For DC, we adopt DC2 
because of the success in explaining the trend of Ry in terms of 
the UV slope and the 2175 A bump (Section r3.2.3t . 

We start from the size distribution of grains processed in DC2 
for 10 Myr. Then, we apply the condition of WIM. In Fig. [5] we 
show the results at t — 3 Myr and 5 Myr in WIM. Around 5 Myr, 
the number of small grains is recovered to the level of the MRN 
distribution. In other words, if grains pass their lifetimes in WIM 
more than in DC, the grains are shattered too much to be consistent 
with the MRN distribution. This implies a short lifetime of WIM. 
Combining this short lifetimes of WIM with a theoretically implied 
timescale for the phase exchange (a few x 10 7 -10 8 yr; Section |231 l, 
we obtain a picture that a large fraction of warm medium is in a 
neutral form and a certain small fraction is ionized. It is interest- 
ing to point out that such a short timescale is consistent with the 
recombination timescale as mentioned in Section [3.2.1l 

The corresponding extinction curves are shown in Fig. [6] The 
Milky Way extinction curve is indeed recovered by the phase ex- 
change. This demonstrates that it is really possible to reproduce the 
Milky Way extinction curve by considering dust grains processed 
in multiphase medium. 

The above phase exchange model is too simple, and the realis- 
tic ISM has more continuous density distribution and m ore compli- 
cated structure of turbulence ( Wada & Norman 2001). Such com- 
plexity should tend to eliminate the specific features such as accu- 
mulation of grains around a ~ a few x 10 -6 cm in DC and selective 
grain destruction at a ~ 10~ 6 cm in WIM. Thus, we expect that 
the grain size distribution becomes smoother in realistic ISM than 
we calculate in this paper. 
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Figure 7. The grain size distributions of (a) silicate and (b) graphite after t = 50 Myr (solid line) and t = 100 Myr (dashed line) in WIM for the MRN size 
distribution extending up to a = 10 -4 cm as the initial condition (dotted line). 



4.2 Upper and lower limits of grain size 

According to M RN, the upper grain radius is ~ 0.25 /J,m (see also 
iKim et alj[l994h . Coagulation has negligible influence on grains 
larger than a ~ 0.2 /j,m both for silicate and for graphite because 
they generally obtain larger velocity than the coagulation thresh- 
olds. Thus, if there is no grain with a > 0.2 fim initially, it is not 
possible to make such large grains by coagulation in ISM. 

Even if grains larger than a ~ 0.2 /j,m form by condensa- 
tion in stellar ejecta , shattering could destroy such large grains. 
iNozawa et al.l J2003h show that silicon grains with o > 0.2 ;um 
form in Type II supernovae. In the outflows from evolved late-type 
stars, the grain radius i s expected to become of order ~ 0.1 /im 
( iGail & Sedlmav3l 19991) , and grains with a > 0.2 /im may have a 
chance to form. It is interesting to note that grains with a > 0.2- 
0.3 fim are accelerated above the shattering threshold in WNM. 
Thus, shattering in WNM may play a central role in determining 
the upper limit of the grain size in ISM. 

In order to examine whether or not shattering in WNM really 
plays a role in determining the upper limit of the grain size, we per- 
form a test by adopting an initial grain size distribution extending 
up to a — 1 /im with the total mass of grains conserved. Then the 
evolution of the grain size distribution is calculated by applying the 
conditions in WNM. Fig. [7J shows the results. We observe that the 
grains with a > 0.2-0.3 fim are significantly shattered in 50 Myr. 
Thus, shattering in WNM is a strong candidate for the determining 
mechanism of the upper limit of grain size. 

In Fig. [8] we show the corresponding extinction curves. The 
initial extinction curve is significantly lower than the observed one 
because large grains tend to have low mass absorption coefficients. 
However, after 50 Myr, the level of the extinction is already consis- 
tent with the Milky Way curve. This means that shattering of large 
grains in WNM is efficient enough to reproduce the upper grain 
size consistent with the observed Milky Way extinction curve. 
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Figure 8. The extinction curves calculated for the size distributions in Fig. 
UJ The dotted line represents the initial extinction curve (MRN size distri- 
bution extending up to a = 10 -4 cm). The solid and dashed lines show 
the extinction curves of grains at t = 50 Myr and t = 100 Myr in WNM, 
respectively. 

4.3 In the context of galaxy evolution 

The efficiencies of shattering and coagulation are affected by the 
grain abundance. This indicates th at metal-poor ga laxies, which 
are generally poor in dust content dlssa etal.| [l990). have differ- 
ent grain size distributions. Here we examine the metallicity de- 
pendence of shattering and coagulation. We assume that the dust- 
to-gas ratio is proportional to the metallicity Z\ that is, we adopt 
TZ = 4.0 x 1O" 3 Z/Z and 3.4 x 10" 3 Z/Z Q for silicate and 
graphite, respectively, in equation J2 1 1 . In other words, the dust 
density in the ISM is proportional to the metallicity, and we expect 
that the effects of shattering and coagulation become weak as the 
metallicity decreases. The turbulence model and the grain veloci- 
ties are not changed, which means that we implicitly assume that 
the parameters listed in Table[T]are fixed. 
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Figure 9. The grain size distributions of (a) silicate and (b) graphite at t = 10 Myr in WIM. The solid, dashed, dot-dashed lines show Z = 1, 1/3 and 1/10 Zq , 
respectively. The dust abundance (vertical axis) of the dashed and dot-dashed lines are multiplied by 3 and 10, respectively, to offset the low dust abundances. 




Figure 10. Same as Fig.|9]but in DC2. 



We test WIM and DC, where shattering and coagulation, re- 
spectively, are the most efficient among the various phases. In 
Fig-El we show the grain size distributions in WIM at t — 10 Myr. 
We apply a longer timescale than adopted in the other part of this 
paper to enhance the effect of shattering. We observe that the shat- 
tering effect is significantly reduced at 1/10 Zq. The same is true 
for coagulation in DC as shown in Fig. [10] where we adopt DC2 
because of the success in reproducing the trend of Rv in terms of 
the UV slope and the 2175 A bump (Section [3.2.3l l, Thus, as the 
metallicity decreases, the relative importance of processing by in- 
terstellar turbulence becomes minor in determining the grain size 
distribution. This indicates that the initial grain size distribution at 
the grain formation in stellar ejecta is relatively preserved in metal- 
poor galaxies (typically Z < 1/10 Zq), although we should keep 
in mind that there are other processes, such as interstellar shocks by 



supernovae, which could modify the grain size distribution in any 
metallicity. 



4.4 Toward the grain evolution in protoplanetary discs 

The condition of turbulence in the circumstellar discs is still 
unclear. Let us con s ider protoplanetary discs. According to 
iNomura & Nakagawal d2006l) . turbulence is very weak (SV ~ 
0.01-0. lc s , where c s is the sound speed). The acceleration by tur- 
bulence will be marginal in this case, and grain motions are more 
likely to be Brownian. As a result coagulation is at least as efficient 
as in DC. As shown in Figs.|2]and[3] small grains with a < 10 -6 
cm are strongly depleted in DC because of coagulation. Thus, we 
can justify that the grain size distribution in protoplanetary discs 
is biased to radii > 10~ 6 cm. Moreover, because grain veloci- 
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ties are expected to be lower than the coagulation threshold even 
at a > 10 -6 cm, gra ins grow fu r ther. 

As shown by ISano et ID (2000), the grain size in proto- 
planetary discs is important in determining the unstable regions 
for magnetorotational instability, which induces MHD turbulence 
tealbus & Hawlevll 19981) . Consequently the grain size distribution 
is further affected by the presence/absence of the turbulent motion 
determined by the instability/stability condition. The coupling be- 
tween turbulence and grain size is interesting to investigate as a 
future work. 



5 SUMMARY 

We have investigated the effects of shattering and coagulation on 
the dust size distribution in turbulent ISM, adopting the typical ve- 
locities of dust grains as a function of grain size from YLD04. By 
using a scheme of grain shattering and coag ulation which we hav e 
developed in this paper based on JTH96 and lChokshi etafl Jl993h . 
we have calculated the evolution of grain size distribution in turbu- 
lent ISM. Since large grains tend to have large velocities because 
of decoupling from small-scale turbulent motions, large grains tend 
to be shattered. On the other hand, because of small surface-to- 
volume ratio, large grains require more time to be destroyed. 

Large shattering effects are indeed seen in WIM for grains 
with a > a few x 10~ 6 cm. In the supernova shocks, such small 
grains are decelerated quickly by gas drag and larger grains tend 
to be shattered more efficiently (JTH96). Graphite grains are pre- 
dicted to be shattered also in CNM, but the result in CNM is sen- 
sitive to the threshold velocity for shattering. Coagulation signifi- 
cantly modifies the grain size distribution in DC. In fact, the cor- 
relation among Rv, the carbon bump, and the UV slope in the ob- 
served Milky Way extinction curves is qualitatively reproduced by 
the coagulation in DC. We have also shown that the upper limit of 
the grain size in ISM can be determined by the shattering in WNM. 

If a large fraction of ISM experiences either WIM or DC, the 
grain size distribution in ISM may be determined by a balance be- 
tween shattering in WIM and coagulation in DC. Considering that 
the effects of shattering and coagulation become small in metal- 
poor environments, the regulation mechanism of grain size distribu- 
tion is quantitatively different between metal-poor and metal-rich 
environments. 
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